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Abstract 

We study the hadronic form factors of r lepton decays r — > Kir(n)v. We compute one loop 

i 1 

corrections to the form factors using the chiral Lagrangian including vector mesons. The countert- 

a? 

erms which subtract the divergence of the one-loop amplitudes are determined by using background 
field method. In the vector form factor, K* resonance behavior is reproduced because the diagram 
with a vector meson propagator is included. We fit the data of the hadronic invariant mass spec- 
trum measured by Belle by determining some of the counterterms of the Lagrangian. Besides 
the hadronic invariant mass spectrum, the forward-backward asymmetry is predicted. We also 
study the effect of CP violation of a two Higgs doublet model. In the model, CP violation of the 
neutral Higgs sector generates the mixing of CP even Higgs and CP odd Higgs. We show how 
the mixing leads to the direct CP violation of the r decays and predict the CP violation of the 
forward-backward asymmetry. 
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I. INTRODUCTION 



The tau hadronic decays are unique as t 
CP violation beyond the standard model 



decays has been searched in r — > -k-kv 13 1 



re decays can be useful for the search for the new 
lNl2l] . CP violation in tau lepton semileptonic 
and t — > Kiru 



[14] . flo1 | modes. Recently, Belle 
15[ puts constraint on the CP violation parameter of three Higgs doublet model with their 



latest data. Babar |16l also searched CP violation of r — > K^tt v(n)7i° (n > 0) and obtains 

nn 

non-zero CP asymmetry which sign is opposite to the standard model prediction |17|.|18l]. 
In r lepton system, another CP violating observable, e.d.m. (electric dipole moment) is also 



searched 



19]. 



To predict the direct CP violation of the hadronic r decays, the strong phase shifts are 
important quantities and the quantitative prediction on the strong phase shifts is neces- 
sary when extracting the weak CP violating phases from the experimentally observed CP 
asymmetries (j[ Q, 11 1. This is a reason why we study the hadronic form factors. 

The hadronic form factors for the decay r — > Kiru and SU(2) isospin vector form factors 
have been studied with various methods. The common future of them is the effects of vector 
mesons (p,K*) and higher resonances are included. In Ref. jol. I20]. vector dominance models 
are studied. One loop corrections to the SU(2) vector form factor are studied with the 



resonance chiral theory in 



In Ref. 



22N25| the t — > Ktcv form factors are predicted 



using the chiral theory combined with the dispersion relations. In our previous study, we 
use the resonance chiral Lagrangian including the one loop corrections to the self-energy 



of resonance 



101 ]. We also note in the experimental study 26|, the Breit-Wigner form for 



101 ] . By using the chiral 



several resonances is used to fit the data of the hadronic invariant mass spectrum. 

In this paper, we use different approach from the previous study 
Lagrangian including vector resonance j^], we compute the one loop corrections of pseu- 
doscalar mesons to the form factors. Including the Feynman diagrams with a vector meson 
propagator, one can reproduce the resonance behavior near the poles of the resonances while 
keeping prediction at threshold region consistent with the chiral symmetry. Our Lagrangian 
includes several new counterterms related to vector mesons in addition to the counterterms 
which are present in chiral perturbation theory. The coefficients of the counterterms are 
different from the chiral perturbation theory. 

Since Belle and Babar reported the precise measurements of the branching fractions of 
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t~ — > K s ix~v [26], t~ — > K~tt q v 30| and r~ — )■ K~r\v 3jJ, |32| . we can compare our 
prediction of the hadronic invariant mass distribution with the experimental data. We have 
determined the finite parts of the coefficients of counterterms so that the hadronic invariant 
mass spectrum is reproduced. 

Once the form factors are fixed, one can use them for predictions of various distributions 
within the standard model (SM) and beyond. We first compute the angular distribution 
and the forward-backward asymmetry (FBA) for r — > Kixv and r — > Ki]v in SM [33}] . Fur- 
thermore, CP violation of FBA is predicted with a two Higgs doublet model. In type II 
two Higgs doublet model, within the tree level approximation, the charged Higgs couplings 
with quarks and leptons are written in terms of Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix. However, if we take into account the one loop corrections to the masses of quarks and 
leptons due to the neutral Higgs exchanged diagrams, CP violation of the neutral Higgs 
sector becomes a new source of CP violation of the charged Higgs Yukawa couplings. We 
show the CP violation of the type II two Higgs doublet model can be probed by the direct 
CP violation of the r hadronic decays. 

The paper is organized as follows. In section II, we show the hadronic chiral Lagrangian 
with the vector resonances which are relevant for the form factors. The counterterms are 
also given. In section III, we compute the form factors. In section IV, we calculate the 
hadronic invariant mass spectra of the decays r — > Kt^{j])v. The spectra are compared 
with the experimental data and the FBAs are predicted. In section V, we explain how CP 
violation in neutral Higgs sector reveals itself in the charged Higgs Yukawa couplings in a two 
Higgs doublet model. We also calculate the CP violation of FBA of the hadronic r decays 
and the numerical result is presented. Section VI is devoted to discussion and summary. In 
appendix, we give some details of the derivation of the formulae used in the text. 

II. CHIRAL LAGRANGIAN WITH VECTOR MESONS 

The leading order 0(p 2 ) of chiral Lagrangian with r] — r] s mixing term and vector meson 
mass term is given by, 

f 2 

C = ^-Ti(D L UD L U j ) + BTr{M(U + U j )}-ig 2p Tr^MZ-^M^)-rio 



1 M 2 
+ -<9^ <9% - -fvl + MyTr 



(1) 



where U is the chiral field which is given as U — exp(2iir/f) = £ 2 . tt is SU(3) octet 
pseudoscalar and B is a constant parameter. r] is U(1)a pseudoscalar of which mass is 
denoted by M and g^ v is the coupling for r] — r] 8 mixing. The covariant derivative for the 
chiral field U is given by, 

D Lfl U=(d lt + iA L jU, (2) 

where the external gauge field of SU(3)l denoted by Al is introduced. is the vector 
nonets and a M is defined as, 



a, 



2i 



(3) 



The form of the mass term of vector mesons is identical to the that of the unitary gauge 
fixed version of hidden local symmetry approach 22-2^|. The kinetic term of the vector 
mesons is not included in the leading order. This treatment is important when including 
loop corrections in a systematic way. Note that we have added the chiral breaking term 
by M = diag(m u , m^, m s ) for the pseudoscalar mesons. The chiral breaking term \ i n the 
isospin limit can be written in terms of the masses of tt and K mesons as, 

ABM 



X 



P 



m~ 
ml 









2m\ 



mz 



(4) 



T) — T) mixing 



When computing the form factors for r — > Kr/^v, they are sensitive to the mixing angle 
of r) and t]', We first summarize the mixing of the octet and singlet pseudoscalar meson at 
one-loop order. The self-energy correction for rjo and rjs sector in one loop is computed with 
the interaction terms shown in Appendix D, 



V8 VO 



— 1 2 



m: 



SMl 8 M n 2 s + 6M& 



08 



M 2 8 + 5M 2 



08 




(5) 
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where the tree level mass squared matrix elements M 88 and Mq 8 are given by 4my 3 m "* and 
^^{m 2 K — ml) respectively. z 8 8 — 1> ^^08' anc ^ ^^fs are one 1°°P corrections and they are 
given by, 

z 88 = 1 + 6c/i A - - — ±(2m 2 K + ml) - 



8ml- — 5m 2 , 2 2 2 
+ 2(4m x - 3mJfx K + 6/j, n m w 



08 " ^ L 

+ 4(2m 2 K + m 2 n )t r 3 + 8m 2 K t r 5 

2 16m 2 K - 7m 2 n 16m 2 K - 6ml 2 
5M 88 = g Vvs + 3 m + 

2 92 P f * ■ a ( ,5m 2 -8m| 

H 5— COS 6>08 Sin 6*08 (/V ~ A^j 



B ~^ uo u — " ™ 3V3 

iaTr ( 2m K + ™l)(4m 2 K -ml) mt + 2(2m 2 K - ml) 2 
+ 16L 6 3^ + 32L 8 — , (6) 

M 



where c = 1 — -^p and /2 P denotes 327r ^ 2 l°g ~- f> (-P = ^1 K, T], 77'). /i denotes renormalization 

scale. We also introduce the notation; fi Va = /i^cos 2 ^ + sin 2 6*08- ^(z = 3,5) and 

L\{i = 4,5,6,8) are the finite counterterms which are defined in Eq. fl2T|) and Eq. fl23|) . 9 08 

denotes the mixing angle at the leading order and is given by, 

1 2\M 2 I 

9 08 = -2 arctan Mj3%- ( 7 ) 

The self energy in Eq. (j5"]) can be diagonalized with the following transformation, 





H 


' y/I^ 1 


l( 






V 11 





cos ^08 sin ^os 
- sin ^os cos 6*08 




(8) 



where, 



(9) 



m 2 
m 2 , 

cos6» 08 -sin0 O8 \ / Mf 8 z 88 + <5M 8 2 8 M$ 8y /z^ + <5M 2 8 \ / cos6» 08 sin6» 08 

sin ^ 8 cos^8 J \M$ 8 y/z^ + 5M$ 8 M 2 J ^-sin0 Q8 cosmos 

We use the transformation Eq. ([8~j) when we compute the form factors for r — > Kr^v decays. 
From Eq.Q, M and M 08 are written by, 

Ml = m 2 + m\, - M 2 88 z 88 - 5M 8 2 8 , 



zV 8 M 2 08 + 5M 2 08 = -^M 2 m {M 2 88 z 88 + 5M 2 88 ) - m 2 m 2 ,. (10) 
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Eq.( !TO|) can be used to obtain the input values for Mo and Mo$ from the mass spectrum 
when the finite counterterms are given. The mixing angle #os including the correction is also 
given by, 



flos = -^ arctan gp» XT 1 "7% • (") 



1 2|M 2 8v ^8i + ^M 2 

2 arCtan M 2 - Mf 8 z 88 - 5M| 



When we compute the form factors related to rj and rf in one-loop order, the mixing angle 
#08 implies one-loop corrected one. The treatment is consistent with the rigorous one-loop 
computation and the difference is at two loop order. 

B. vector meson sector 

Now we turn to the vector meson sector of the Lagrangian. The quantum corrections to 
the chiral Lagrangian with vector mesons have been discussed in several works 



21 



28 



29). 



They study the vector meson loop correction in addition to the loop correction due to 
pseudoscalar mesons. Our aim is to construct the effective theory which can be used to study 
the process in which a single vector meson can be nearly on-shell. The corresponding energy 
region for hadronic invariant mass is Eh < 2My. For the hadronic r decays, this approach is 
valid in the energy region Eh < 1400 ~ 1600(MeV). In the region, the vector meson does not 
contribute to the loop diagram and only the light pseudoscalar mesons loop should be taken 
into account. The loops of the soft pseudoscalar mesons can be systematically included using 
the momentum and loop expansion. We regard the vector meson as the classical background 
field. As for pseudoscalar mesons, we split them into hard classical background field and the 
soft quantum fluctuation. For example, the decay products Kir of K* meson decay have the 
hard momentum Eh ~ My and they are treated as the classical background field. Though 
the vector mesons do not contribute in the loop diagrams, they contribute to the amplitude 
as intermediate dressed propagator which connects the 1 PI (Particle Irreducible) vertices 
of vector mesons and pseudoscalar mesons. The self-energy of the vector meson and 1 PI 
vertices with or without vector meson legs can be systematically improved by taking the 
quantum corrections of the soft pseudoscalar meson loops into account. 1 PI vertices can 
be renormalized by adding the counterterms. What kind of counterterms should be added 
depends on the number of the pseudoscalar loops N and the number of the external legs of 
the vector meson Ny. We focus on the chiral limit. The number of the derivatives for the 
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counterterms is determined by superficial degree of the divergence of the 1 PI diagrams. As 
we prove later, the superficial degree of divergence of 1 PI diagram of N loop order and Ny 
external vector meson legs is given as, 

u = 2N + 2 - N v . (12) 

This formula tells us the types of the counterterms which should be added when we carry out 
N loop order computation. In general, the local counterterms and the finite counterterms 
can be classified with the number of derivative and the number of the vector mesons ny 
in the Lagrangian. The interaction term with n d derivatives and with ny vector meson fields 
has the form of, 

£ = F(i)d nd V nv , (13) 

where the Lorents indices are contracted appropriately. F(£) denotes some function of the 
chiral field. The derivatives can act on both the chiral field and the vector field V. Since 
the number of derivatives of the vertex of counterterms is equal to the superficial degree of 
divergence u, the divergence of the N loop order Feynman diagram with Ny external vector 
mesons can be subtracted by the counterterm with the following number of the derivatives 
and the vector meson legs, 

(n d ,n v ) = (2N + 2-N v ,Ny). (14) 

In table I, we show (nd,ny) for a given set of iV and Ny. The lowest order Lagrangian 
corresponds to iV = case in the table I and it includes the interaction terms of the type, 

(n d ,n y ) = (2,0),(l,l), (15) 

where (n d ,n v ) = (1,1) corresponds to the term of Tr[V M a: M ] in Eq.([T]). The lowest order 
Lagrangian includes mass term of the vector mesons while the kinetic term is not included. 



This is in contrast to the approach of 



28 



291 ] where the vector boson is treated as gauge 



boson and the kinetic term is included in the leading order Lagrangian. The vertices of only 
vector mesons do not contribute to vertices in any 1 PI loop diagrams since the vector meson 
does not contribute in the loop diagram. Such vertices include the tree level mass term with 
(nd,n v ) = (0,2) and N loop order counterterms with {n^ny) = (0,2iV + 2) which has the 
form V 2N+2 . Now we prove Eq. (1121) . We consider a 1 PI diagram of N loop order. N loop 



N 
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2 


N 


nd + ny 


2 


4 


6 


27V + 2 


Ny 













(2,0) 


(4,0) 


(6,0) 


(27V + 2,0) 


1 


(1,1) 


(3,1) 


(5,1) 


(27V + 1,1) 


2 


(0,2) 


(2,2) 


(4,2) 


(27V, 2) 


3 




(1,3) 


(3,3) 


(2JV-1,3) 


4 




(0,4) 


(2,4) 


(27V - 2,4) 


5 






(1,5) 


(2JV-3,5) 


6 






(0,6) 


(27V -4, 6) 


27V + 1 








(1,27V + 1) 


27V + 2 








(0, 27V + 2) 



TABLE I: The number of the derivatives and the number of the vector meson external legs ny. 
(rid, n v) a r e shown for a given set of numbers of loop TV and the number of the vector meson legs 
N v . 

order diagram includes the N loop diagrams with the tree vertex as well as the diagrams 
with higher loop order vertices with the number of loop n L less than TV. We denote the 
number of n loop order type vertices with rid derivatives and riy vector meson fields included 
in the diagram as n^[rid, ny}- Note that rid + ny = 2n + 2. The total number of the n loop 
order interaction vertices in the 1 PI diagram is given by 

2n v + l 

Ni n) = Yl n^[2n + 2-n v ,n v }. (16) 

n v =0 

Although in N loop order 1 PI diagram consists of the various loop order vertices, the number 
of the vertices must satisfy the following relation 

N 

N = n L + ^nNi n \ (17) 

71=1 

The number of pseudoscalar meson internal propagator Ib is written as, 

N 

I B = Y,Ni n) +n L -l. (18) 

n=0 
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Then one can compute the superficial divergence u of the 1 PI diagram, 

N 2n+l 

u = 4n L -2I B + J2 ^(2n + 2-n y )n(, n) [2n + 2-ny,n y ]. (19) 

n=0 ny=0 



The last term of Eq. ffl9|) is the number of the derivatives of the diagram. Substituting 
Eq.fflgjl and Eq.flU]) into Eq.(p3]), one can show Eq.(JT2P as, 

JV 2n+l 

u = 2N + n v n{ v ] [ 2n + 2 - n v ,n v ], 

n=0 ny=l 

= 2N + 2-N v . (20) 

In the same way as the chiral perturbation theory, we rely on the momentum expansion. 
Because the loop momentum of the pseudoscalar mesons is soft, the expansion is valid. In 
the Lagrangian, there is no kinetic term for the vector meson at the leading order. The 
kinetic term is generated as the loop correction of the pseudoscalar mesons. According to 
Eq.( fl2|) . one can extend the chiral counting to the case with vector mesons. In generalized 
chiral counting, the vector meson field is counted as 0(p) and the chiral breaking x is 
counted as 0(p 2 ). The couterterms for 0(p A ) are obtained by computing divergent part of 
the one loop corrections due to pseudoscalar mesons. We use the background field method 



and the corrections can be computec 
are consisitent with chiral symmetry 



and the counterterms can be determined so that they 



34| . The outline of the derivation is shown in Appendix 
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A and they are given by, 

L c = K 1 2iTr(a± IJi a± u )(D l j i v v - D u v^ + i[v^v v \) 

- ~ {K 2 Ti{^F L ^){D^ u - D v Vp + v v \) + K 3 Tt(D^v u - D v v^ + i[v^, v v ]) 2 ) 

+ (fQTr{(£M£ + gMg)v 2 } + K 5 Tr{M(U + rf)}Tr(y 2 )) 
+ K 6 Tr(v p a^)Tr(v p axn) + K 7 Tr(v 2 a ± ^a p L ) + K 8 Tr(v 2 )Tr(a ± ^a^) 
+ K 9 {Tr(v 2 )} 2 + K 10 Tr(v 4 ) 

+ irtfrMTtiiZMZ - eWK} + ivoTMTr{M(U - rf)}Tr(v 2 ) 



+ 



T 3 ^i^ Vo TrM(U + rf)TrM(U - f/ f ) + T 4 f^j Vo 2 (TtM(U - 



+ iT 5 ^^ri Tr(MUMU - MU^MU^) 
f f 

+ T 6 rj 2 Tr(MUMU + MU ] ' MU ] - 2M 2 ) 

+ i^ Vo [T 7 Tr{M(D Lfi U(D^ L UyU - U^D L ,U{D^)} 
+ T 8 Ti(M(U - rf))Tr(D L rU(D£U)^ 



+ L^Tr^C/p^)} 2 + L 2 Tr{D£U(D u L Uy}Tr{D Lfl U(D Lv Uy} 
+ L 3 Tr{D£U(D L ,UyDlU(D Lu U) j } 

+ L 4 Tr(D LlM U(D£uy)Tr{M(U + C/t)}^ + L 5 Tr{D Lll U(D^U)^(UM + MU^)}^ 



+ —{L 6 {^M(U + C/t)} 2 + L 7 {TtM(U - C/t)} 2 } 

1 Q R2 

+ L 8 ——Tr(MUMU + MU ] ' MU ] ) + iL 9 Tt{F l ^D"U(D u U) 1 } 

+ HiTrF^Fj? + H 2 Tr(M 2 ) + H 3 M 4 , (21) 

where u M = ^p(V^ — The covariant derivative is defined as; D^v u = d^v u + i[a^,v u }. 
a± in Eq. (j2"Tj) is given as, 

oa„ - . (22) 
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fcl = 1 


h = -6 


r = 2c 2 +l 

OZ 


Ai = -i 

o 


k 2 = 1 


t 2 = -2 


r _ l+2c 2 
1 2 — ifi 


r,2 f2 

A _ 5 | 92p/ 


k 3 = 1 


J 3 — — Ts 


10 


A 3 = 1 

z 


z 


+ — 11 
t-4 — o" 

y 


^4 = | 

o 




fc 5 = \ 
z 


*5 = — 1 



T 3c 

1 5 — "a" 

o 




fc 6 = 4c 


^6 = — I 

o 


„2 f2 
r 11 9-2pJ 

-^6 1/1/1 Q/1 D!2 




ky = 6c 


j. 3c 
17 — — y 


r 7 = o 




fc 8 = 2c 


tg = -f 


„2 j-2 

■p 5 | 92 P J 
L 8 48 "i" 




fcg = -3 




r 9 = 4 




fclO = -3 









TABLE II: The coefficients of the counterterms; fc^Tj and Aj. c = 1 yw. 



The coeffcients of the counterterms are splitted into the finite parts and divergent parts as, 

Ki = \k i + K r t {i = l ~ 10), 
Ti = \U + T[(i = l~6), 
U = XT i + L r i (i = l~9), 

Hi = AA i + H < r (i = l~3) > (23) 

where, 

X = -^(C uv + l-\og^), (24) 

with C^y = - — 7 + log47r. The coefficients fc^t^Tj and Aj are given in the Table II. 
From Eq.( l2TT) . we extract the effective counterterms which are relevant for the calculation 
of the form factor of r — > Kttu decay. The effective counterterms which subtract the 
divergence of the amplitudes which contains a vector meson in Fig.l, can be deduced from 
the counterterms shown in Eq. (|2"T|) . They are the counterterms for the self energy of vector 
mesons, V — > PP vertex, and the production amplitude of the vector meson; — > V and 



11 



are defined as, 



1 



2 — v — \- Vfiu 



+ i^Tr(F^d„nd v n)+C 4 Tr(F^Fl v ), 



C 2 Tr (M±i«l) Tr 



a 



M\2 



(25) 



where F Vflu = d^V v - d u V^ and F^ = d^A Lu - d u A Lfl .. Z v and Q (i = 1, ... ,4) are 
renormalization constants. The coefficients C{ can be written in terms of the coefficients of 
the counterterms of Eq. ( 12T1) , 

r2 \ 2 



d = 2K A 

C 2 = 2K 5 
M 

c 4 



M 



2 \ 2 



2 \ 2 



M 2 /If 2 



'2<r7 2 ' 



/If 2 /If 2 



49/ 2 ' V/ 2 
The finite parts of the counterterms also satisfy the relations similar to Eg. (126 



(26) 



C{ = 2K\ 



Ml 



v 



2gp 



CI 



c: 



2K\ 

My 



2 \ 2 
V 



gf 



(K\ - K. 



. M? 



v 



3 V/ 2^ 



/If 2 /If 2 
\ K 2 - K 3] 



(27) 



Agf 2 V ' ^ 2/2 ,- 

One can extract the counterterms for the 1PI vertex of the type — > PP. They are given 
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f 

+ %{-K A 



1 / M 



2 \ 2 



2/ 2 WPJ ' ~f 

2 \ 2 



+ 4^}(TrA L ,{[n,d^], X }) 



L 4 K 5 ( My 



P P \2g 2 P 
+ ¥*TrA Lll {&>K, [tt, ; 

where C5 and its finite part C5 are given by, 



c 



M v ( K (r) _ jXr) M v (r) (r) 



(r) 



(28) 



(29) 



III. THE FORM FACTORS AT 0(p 4 





FIG. 1: Feynman diagrams which contribute to the form factors. The diagrams are classified into 
two categories , one of them corresponds to 1 particle irreducible diagram. The other is the diagram 
with a vector meson propagator. The crossed circle denotes the weak vertex. The circle denotes 
the interaction vertex of the vector meson and the two pseudo scalar mesons. 



In this section, we compute the form factors. The matrix elements of the current u^/^s 
are obtained by identifying the quark current with the corresponding current of the chiral 
Lagrangian, Eq.((T]) and Eq.( l2~T]) . In r — > Kttu process, K* meson which is produced by 
the strangeness changing current can contribute to the vector form factor. The resonant 
contribution is significant when Kir invariant mass is near to the resonance pole. In our 
framework, the resonant contribution is included in the second diagram of Fig.l. We take 
into account the resonance contribution by using the vector meson propagator with one-loop 
corrections to self energy. Since the propagator have a pole in complex plane, the effect of 
the width of resonance is also included. Thus we can reproduce the resonance behavior. 
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There are three parts of the diagrams of Fig.l. The first one is 1 particle irreducible (1 PI) 
diagrams and the diagrams with one loop corrections are shown in Fig.2. They correspond 
to the one loop corrections to W + — » K + tt° vertex. They include all the diagrams which 
are also present in chiral perturbation within one loop. Their contributions to the matrix 
element (K + TT°\ui'j fl SL\0) become, 



where Q = Pk + Pn and q = Pk — Pn- ^pq denotes the mass squared difference 

Apq = mp — trig. Spq denotes the sum of the mass squared £pq = mp + itlq. The 

loop functions Ip,x® P , J9 P are given in Eq.f lBip and Eq. flB3j) . and z n are finite wave 
function renormalization and they are given as, 



(K + n \uZ^s L \0) 



ipi 




+counterterms, 



(30) 



z K = 1 + 3c(}x K + -(n m + - —{L\{2m 2 K + m%) + L r 5 m 2 K }, 

g 

z n = 1 + c{2fi K + 4^) - — {L\{2m\ + m\) + L\m\}. 



(31) 
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Including the finite part of the counter terms, the result of the 1 PI part is 
(K + 7r°\u Lltl s L \0) 



ipi = 271^ -^{H K K + H Km ) + ^^{lUii K + 3/^ + 11/^) 



3 f M< 



v x 2 
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2 ' ~ T7 2 

2m\ + m\ M v 



{Hrtt + H Kr}8 + 
r2 ™2 



+ 



2g*p\2gip 4 p 
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2g 2 f 



2V2Q 2 



'1 - 
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2 f 2 ' 



A '--% {5c g2_ (5c _ 3)SK7r} 



^■KqJKr] COS 2 # 8 + Ak v 'Jkt]' sin 2 # 8 

2 A2 

1 JKtt + 

s s 



+ ^(1 



{3cQ 2 -(3c-l)S^}} 
A 2 



+ 



]f 2 V 2^2/ 2 
1 AxttQa* 



2V2 Q 2 

-2 \ 2 



My 



c 2 3/i m + 2/2 K - 5/^ 
4^ + V/ 2 



-Jx^sin ^08} 



(10/iK + 3/i^ + 11^) 



3 / M 2 V Q 2 



M 2 



2<? 2 / 2 



x 



r E K77 + 2m 2 K + ml M 2 _ g 



2^ 2 / 2 4 / 2 5 P P 2g 2 P 5 

and C|" are finite parts of the counterterms L, and C, respectively. The function Hpq is 
written in terms of the functions defined in Eq.( lC9j) as, 

1 



(32) 



Hpq = j^(Q 2 M PQ - L PQ ). 



(33) 



Jpq can be found in Eq. flClOp and Eq.( 1C12l) . We also introduce the following notations in 
this paper. 



Y Km = Y Kri cos 2 9 0S + Y Kn > sin 2 6> 08 , 



where Y = H, M r and L which also appear in the following equations. The diagram with a 
vector meson propagator is shown in Fig.l. It includes the diagram with a K* propagator, 
K* K + tt° vertex and W + — > K* production amplitude. The self-energy of the propagator, 
the vertex and the production amplitudes include one-loop corrections. We first compute 
the one loop corrections to K* — > Kit vertex which are shown in Fig. 3. 
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FIG. 2: The one loop 1 particle irreducible Feynman diagrams contributing W 4, 
factors. The counterterms are shown by the black blob vertex. 



K + ir° form 



\AA. 






FIG. 3: The one loop Feynman diagrams contributing to K* + — > K + ir° vertex. The black blob 
denotes the counterterm. 
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2 



+ -^-p%\S zrz-k — 1) + counterterms 
i / - i /-' 

lf^T^lT/* 2 + 0(10^ + 3^ + 11/1,) 

- 32L ; 2 < + m --16L- E ^ 



/2 5 /2 

+ -^pi C H 2m K + ml) + CJmy + ||(QV - A^g,) 

+ g J 4 ~~{^'A'7r(AA: 7r JjiTTr — -A Km J K , ns ) 

2 - 3 - 

+ c(Q — S A - 7r )(A^ 7r — -A^ g Jftr^g)} 
= {-2^(! - MsW 3 *- ~ Z v) ~ wY Q ^ ~ Q%) 



2g2f 2^ V, g/2 



r M 2 i 

9 y/ , {(^. + g 2 5i?^)(l - 2^3) - - C^m 2 , + m 2 )} 

M 2 



+ IQgp ^y-^K + Z^tt + /%J + C(10//K + 3/i^ 8 + ll/vj 



. r 2m^ + m 2 



- 32L" — ^- — £ - 16L 



4 /2 5 /2 

My Q - 1 - 2 

+ g~j4 ~{^JC7r(Ajir7r-7g7r ~ "(A/^ Jft^ COS #08 + A^/ Jjc^' Sin 6*08)) 

3 

+ c(Q — S A : 7r )(A^ 7r J^ 7r — -(A Kri J Kri cos O 8 + A^/ J^/ sin 6>os))}, 

(34) 

where A Kr]s J Km = A Kv J Kri cos 2 9 08 + A Kr) i J Krj i sin 2 # 08 . Next, the propagator of the if* 
meson is obtained by including one loop self energy corrections. Using Eq.( 1C2j) . the K* 
meson propagator is given by iD pp where D pp is given by, 

_ Q„Q P SB 
_ M^+SA+QHB 

where the self energy corrections 5 A and SB in this section are identical to the K* mesons 
ones given in Eq.( 1C8j) . 

5 A = 5A K *, 5B = 5B K *. (36) 
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The K* production amplitude with the one loop corrections are shown in Fig.4 and is given 



i i 
i i 



\/\/\ 



FIG. 4: The one loop Feynman diagrams contributing W + — > K* production amplitude. 

by, 



^[C im 2 K + C 2 {2m 2 K + ml)]g^ + V2C 4 {Q P Q U - Q 2 g pv ) 



V2g l 
3M 2 

i r r w 2 i 

-j=- ^5 A + Q 2 5B - -^jf(L K7T + L Kim )j 9fU , 

+ |^ - 2gC[ -5B + -^f (M r Kw + M r K J j (QV - Q M Q„) 



where, 



(37) 



C * = I 1 " ^ ' + 1 - } + ^ (38) 



My ( _ Ml 



J®f is defined as, 



+ ^ (^f^ - 2AEpq - 2(„ e + . (39) 

is defined as; JJH K = Jff cos 2 Oqs+J^ sin 2 # 8- Now one can assemble the contribution 
from the diagram with a K* propagator to the form factor. One can write K* production 
amplitude of the weak vertex and K* — > Kn decay amplitudes as, 

(K;\uZ^s L \0) = g v ,G + (Q 2 g Ufl - Q U Q P )H, 
T p (K* + K + n°) = Eq p + FQ p A Kn , 

(40) 
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where G,7i, ,E and T are given as, 

1 1M 2 
G = ^~ g {M 2 v + SA+ Q 2 5B - -^(L K7T + L Km )}, (41) 



1 

n = ^{Zy - 2gCl - 5B + (M^ + M r K J}, ( 1.2) 



M 2 2t7i 2 -\- m 2 Sr' 

+ —^{-3(2^ + fi w + fi m ) + c(10fjL K + 3fx m + 11^) - 32L^ y * - 16L£^f} 



+ {^a-^x^-^ + ^w 2 , (43) 



'2M 2 ^ 2^ 2 / 2 



M£ 1 - Jk v A Kv cos 2 6»os + Jkv'Akt,' sin 2 6> 08 , 

^2 ^ \/t 3(Jk t; A^cos # 8 + ^k??' A Kri < sin 6> s) x UA s 
+ c(C^ - L, K n){J K n — Jj. (44J 

4^Kn 

Using the form factors, we obtain, 

(K + ir°\uEl»s L \0) |^ = t(Eq p + FQ p A K7T )iD^ (g^G + (Q 2 g afl - Q a Q p )H) 

~ My + 5 A [Qfl Q 2 W g 2 ^ M 2 + <L4 + Q 2 5B ' 1 J 



The vector form factor and the scalar form factors are defined as, 

(K + 7r°\uZ^s L \0) = ^{F v (q, - ^Q,) + F S Q,}. (46) 
Then the contribution to the form factors is given as, 



b s 



T? K * 



(47) 
(48) 
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where, 



-iipj 
v 



My 

2g 2 P 



) + 



+ 



V2 



3c (H + 77 ^ cM v 



{10/lK + 3/ir, 8 + llflir) 



H m ^(h +h + 



2/iA' + /Ar + A*r/8 



2 / 2 



M 2 



K 4 — -4L 5 — 



2g*pl2g*f* 4 / 2 

1 c + g 2 H 



2m| + m 2 ^ M 2 



./' 2 



P 



2g 2 f 



'}]■ 



-2£ 



M 2 



M 2 



2^ M v + 5 A 
1 1 



4* + V2^ 

/ 



(49) 
(50) 
(51) 



(1 " 



My 



2g 2 P' 



Ktt'J Kit 



{5cQ 2 - (5c-3)£^} 



+ 



/ 2 

J Kv cos 2 ^8 + J Kv , sin 2 £ 08 {3cg2 _ ^ _ j 



+ 



3A^ ( 

1 A_ff7T 



8/ 2 

K ^ 2rAL -J^ + ^cos 2 ^Jx n + 



-(1- 



My 



A 2 



My 



sin 2 9 08 J Kv '} 



2g*F' + 4r i A Kw + %g 2 P 



V2 Q 2 

3 / Af 2 \ 2 Q 2 

{ 



(Kfytjc + 3/i 8 + 11/+J 



x 



M V K r m K .rr^Kn 

2g 2 P 4 P 5 P 
A Kn E + Q 2 7 



+ 



2m 2 K + m 2 



P 



M 2 
2g 2 P 
M 2 



2g 2 P 



-2G 



Q 2 M 2 +5A + Q 2 SB 

1 A Kw My 



■G 



My 



A(E + Q 2 F) + V2- (jp 



(52) 
(53) 
(54) 



2V2g Q 2 M 2 +5A + Q 2 SB 
For numerical calculation, we use Eg. (151 ft and Eg. ( 1341 which are obtained by omitting the 
two loop order contribution in the numerators of Eg.( l50p and Eg. (153]) . To compare our form 
factors with the ones obtained by other methods, we show the vector form factors in the 
chiral limit, 



Fh PI =~ 



1 

V2 



(1 



Ml Sl 

K r ){l + 3(1 



2^ 2 / 2 ' 



My 



My 



2V2g 2 P M 2 + 5 A 



2g 2 f 2 
1 + 6(1 - 



Ml } + CIQ 



My 



2g 2 P 



P 

g{Cl-ACl)Q 2 



2M 2 



(55) 
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The self energy correction of vector meson, 5A in the chiral limit can be obtained with 
Eq.tfCl, 



5A 



-Q 2 z\ 



3M V H, 



(56) 



where H is given by taking the chiral limit of 



Q 2 M r 
P 



in Eq.flCH) as, 



H 



Q 2 



i . Q 2 



(57) 



12/ 2L 16tt 2 to /i 2 48vr 2 16tt j ' 

To compare our result with those of the other methods, we examine the case that the vector 
meson dominace (VMD) relation M 2 = 2g 2 f 2 holds. Then the vector form factor is written 
as, 



? iVMD. 
V | chiral limit 



1 

71 



My 



My ~ 3HMy 



ZyQ 2 



{l + g(c;-4Cl)^ 7T }+ - 



2My 



p 



(58) 



The result can be compared with the same limit of the form factor in 

M 2 *e 3Re -(") 
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Q 2 -iM K *T K *{Q 2 



(59) 



The difference of the overall factor — is just due to the the definition of the form fac- 

22] , the chiral loop correction denoted by H 



tors. We observe that in the form factor of 



is exponentiated and appears in the numerator of vector meson propagator while in our 
approach with the vector dominance assumption, the chiral correction appears in the self- 
energy function in the denominator of the vector meson propagator. We also note that the 
finite counterterms generate linear Q 2 dependence in the form factor. They include the 
wave function renormalization constant of the vector meson Zy , and the other coefficients 
of the finite counter terms; CJ, CI and C£. One can also compare our result with that of the 
resonance chiral theory [21| . A difference of the form factor in [2l| from Eq. (158j) of our result 
is that the one loop corrections to their form factor depends quadratically on momentum 
squared Q A . This is due to the second derivatives coupling of the vector meson to two pseu- 
doscalars in their anti-symmetric tensor formulation of vector mesons. In contrast to their 
approach, the vector meson coupling into two pseudoscalar meson coupling includes the first 
derivative. Therefore, the form factor of the present approach depends on Q 2 linearly. They 
also consider the loop contribution of all the resonances while in our approach, the vector 
mesons do not contribute in the loop. 
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IV. NUMERICAL ANALYSIS IN THE SM 



To evaluate the vector and scalar form factors, we fix g, My and the coefficients of the 
counterterms by using the decay constants, masses and widthes of the mesons. We also 
use the hadronic mass spectrum. There are ten parameters, {g, My, Z v , C\, L\, L^}, (i = 
1, • • • , 5) to be fixed. 

From the matrix elements of the axial currents, we obtain the pion and kaon decay 



constants 



fn = f\l- c(2fr + Hk) + 4 ( ml + f l ml Ll + ^Ll)\, (60) 



3c , . „ s . , fml + 2m 2 K m 2 K 



fx = f |1 - j (A* + + fi m ) + 4 ( * ' K L\ + -j^Ll) | . (61) 
Using the ratio of fx/ fm we can write LI as follows, 

^iHf - i -3 (5 ""- 2 ' , *'- 3 '''4 <« 2 ) 

If we assume / = /„., from Eq. (160]) L\ is expressed as, 



L\ = 2 £ 5 (1(2^ + Pk) - *$Ll \ . (63) 
mi + 2m K [A ft 



One can take any renormalization scale p at around K* meson mass. We specifically 



choose the value of the particle data group (PDG) [37j, namely p = 895.47MeV. If 
c(= 1 — M v /(g 2 f 2 )) is obtained, L\ and L\ can be fixed. 

From the imaginary part of the self energy for K* meson in Eq. (lC8p . the decay width of 
K* is given by, 

r r m 2 \ - 1 -MMM ( K \ 2 

r-(AV)- lfaMA „ Mt , (64) 



where is defined in Eq. flClip . Once My is determined, g can be fixed with the decay 



width K* (IV.) and Mx»- The relations among Z v , C[, C r 2 are derived by the conditions for 
the pole masses of K* and p mesons. We define K* and p meson masses as the momentum 
squared (Q 2 ) for which the real parts of the inverse propagators vanish, 

M v + Re[6A K *(Q 2 = M 2 K ,;C\, C r 2 )} = 0, (65) 
M 2 + Re[5A p (Q 2 = M 2 ; C{, C r 2 )} = 0. (66) 
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Solving the above equations, one obtains C[ and C\, 

1 



C\ 



CI 



A 



Kit 



{Z v A K * p - Re[AA K *(Mx*)] + Re[AA p (M 2 p )}} 



2m 2 K + 



where, 

AA K * 

AA r 



3 (Ml 



{Ml - Z V M 2 + Re[AA p (M 2 )] + C[rr£} 



f 2 



4 [gPJ [ Q2(Mk + M ^ s) ~ Lk7T ~ Lkvs + + 2flK + 



Ml 
9ft 



Q 1 ( K + \M r K \ + fl L n + \ m 



(67) 
(68) 

(69) 
(70) 



From the condition for the residue of the vector meson propagator ( 1351) . Z v is written as 
follows, 

dRe[AA K *{Q 2 )] 



Z v 



1 + 



dQ 2 



Q 2 =M%, 



(71) 



We use p meson mass of the PDG value 37J] . For K* meson mass and the decay width, 
we fix them with the hadronic mass spectrum of r — > Ktxv. Instead of using g and Zy as 
the fitting parameters, one can use the decay width Tk* and the mass Mk* ■ 

One can write Zy,C\,C\ and CI in terms of K{, K^-, K£ and L\ with Eqs. (1271) and 
(|29|) . Since K\ is related to Z v with Eq. (127]) and is fixed with Eq. (I7T1) . i£g is already 
determined by Mk* • We note that the form factor of Eq. (I5T]) depends on the combination 



C*3 — 4C4, which is written as, 



4cr 



M 2 



K[ + K r 2 - Kl 



My 



(72) 



9 ft V 1 ' " a ^9 2 f 2 
One also notes that C^ is written in terms of K\ + K^K^ and L\. Therefore we choose 
{T K *, My, M K *, K{ + K2, L r 9 } as fitting parameters in the following analysis 



We fit them by using the differential branching fraction of the experimental data 26 ]. 
The differential branching fraction for KPv{P = tt, 77) is given by, 
dBr (r KPu) 1 G 2 F \V US \ 2 {m 2 T - Q 2 ^ 2 



dJQi 



X 



2 5 7T 3 

2m 2 T 

w 



rrvt 



-Pk 



l)pK 2 \Fr(Q 2 )\ 2 +^\F s Kp (Q 2 )\ 2 



(73) 



where px is the momentum of K in the hadronic center of mass (CM) frame. The differential 
decay distribution for r~ — > K s ii~v is shown in Fig. [5j One can see the peak of K* resonance 
around at JQ 1 ~ 900 MeV. 
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The five parameters are determined by fitting the hadronic mass spectrum in the region 
rriK + < \/~Q 1 ~ < 1665MeV with 90 bins data. We also use the PDG values 37]. 
m f ±, m#o, m,,, m v ',m T as inputs. The set of parameters leading to the smallest x 2 / n -d.f 
value are fixed by 



IV = 48.68MeV, M v = 954.0MeV, M K * = 895.4MeV, 



K{ + K\ 



0.04517, LI 



5.068 x 10" 



(74) 



where the obtained xV n -d.f. is 152.3/85. The other parameters are shown in Table (TTH We 
also note 1 — M v /(2g 2 f 2 ) = 0.2688 for this case. It implies that the relation of the vector 
meson dominance, My = 2g 2 f%, is slightly violated. 



TABLE III: Numerical values of the fitted parameters. 



9 


8.582 


ci 


-0.7772 


LI 


2.265 x 10~ 4 


Zy 


0.8276 


ci - \c\ 


0.1811 


LI 


2.313 x 10~ 3 


C{ 


0.2980 


ci 


-1.516 x 10~ 3 






600 800 1000 1200 1400 1600 1800 
V5" (MeV) 

FIG. 5: The prediction of the decay distribution for r~ — > K s tt~v. The solid line corresponds to 
the prediction of our model. The dotted line and the dashed line correspond to the distribution of 
the vector form factor and that of the scalar form factor, respectively. The closed circles with the 
error bars are experimental data [26]. 
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Table HVl shows the fitted values of Mk*-, ^k* and the prediction of the branching fraction 
of our model. The corresponding experimental values are also shown. We have not shown 
the error for the branching fraction, since the systematic error in each bin is not known. We 
also show the obtained slope parameter for K e3 decay defined in Eq.f lTrJj) . A + is the slope 
parameter given by the linear expansion coefficient of the Kir vector form factor (1470 . 

F^(Q 2 ) ~ F^(0) (l + , (75) 



where, 

ml dF^(Q 2 ) 



Ff*(0) dQ< 



(76) 



TABLE IV: Fitting parameters {Mk* , ^k* ) and the corresponding predictions for the branching 
fraction and the slope parameter. The bottom line denotes the experimental values. 





M K *(MeV) 


rV»(MeV) 


Br(r — >■ Ksir v T ) 


A+ 




894.5 


48.67 


0.4023% 


0.02236 


Exp. 


895.47 ±0.20 ±0.74 


46.2 ±0.6 ±1.2 


(0.404 ± 0.002 ± 0.013)% 


0.02485 ± 0.00163 ± 0.00034 
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FIG. 6: The absolute value of the vector form 
factor. 
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FIG. 7: The argument of the vector form fac- 
tor. We added 180 deg. to argFy. 



We investigate the property of the form factors obtained in the present work. Using the 
fixed parameters, we show the absolute value and argument for the vector form factor in Figs. 
[6] and [71 In the absolute value of the vector form factor |, the effect of K* resonance is 



25 



1.0^ 

0.8 

_ 0.6- 
£ M - 
— 0.4 

0.2 

0.0 



800 1000 1200 1400 1600 
"V^ (MeV) 

FIG. 8: The absolute value of the scalar form 
factor. 
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FIG. 9: The argument of the scalar form fac- 
tor. 



dominant at ^/Q 2 ~ M^*. Furthermore, the effect of K* resonance is seen in the argument 
of vector form factor (argfi 7 ^ 71 "]), because it changes about 180° near \/~Q 1 ~ ~ M^*- This 
property are also seen in Kir scattering j^s)]. However, the behavior of argfF-^] for large 



since our model 



invariant mass region, y Q 2 > 1200MeV is different from the one in 
does not include the higer resonances, i^*(1410) and fC(1790). Figures |8] and [9] show the 
absolute value and argument for the scalar form factor, respectively. The absolute value of 
the scalar form factor is smaller than the absolute value for the vector form factor, since 
there is no K* pole in F£f* as shown in Eq. (I54p . As increasing the invariant mass, the 
argument for the scalar form factor decreases. 

We study t~ — > K~r]i> decay using the parameters fixed with r~ — > K s ir~v decay. The 
form factors for Krj are given in Appendix E. Figure [10] shows the prediction of the decay 
distribution for t~ — > K^rju. It is found that the contribution of vector form factor is 
dominant. The predicted branching fraction for r~ — > K^rju decay is 2.114 x 10 -4 . Since 
the experimental results are Br(r~ — > K^r]v) = (1.52 ± 0.08) x 10 -4 |37l |. our prediction 
is larger than the experimental data. We note that the predicted branching fractions for 
t~ —7- Ks^v and r~ — > K~i]i> decays are 4.030 x 10~ 3 and 1.157 x 10~ 4 respectively with 
the other parameter set of parameters which is obtained by 67 bins data fitting (mx + < 
^fQ 1 < 1400.5MeV). 



We also consider the forward-backward asymmetry 33] for r — > KPv decay. The double 
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FIG. 10: The hadronic invariant mass distribution for r~ — > K~r]i> decay. The solid line and the 
gray solid line correspond to the hadronic invariant mass distribution for t~ — > Kr\v decay and 
that for t — > K s ttv decay, respectively. The dotted line and the dashed line are the vector form 
factor contribution and the scalar form factor contribution of r — > Kr\v decay, respectively. 



differential rate of the unpolarized r decay |6|] is given by 



1 Gl\V, 



F\ v us\ 



[ml 



Q ) Pk 



dyQ^d cos 9 



2 5 7T 3 
+ -A\ F KP\2 



rn 
~4 



mi 




cos 2 9 + sin 2 9 )p 2 K \F« p (Q 



■)2m2 



--p K cos 9Re[F* p (Q 2 )F« p (Q 



V 



-)2\*l 



(77) 



where 9 is the scattering angle of kaon with respect to the incoming r in the hadronic CM 
frame. The forward-backward asymmetry extracts the interference term of the vector form 
factor and the scalar form factor. 

fid cos 9— dBr - ? - » /L;r 

AMQ 2 ) = — 



} 2 dcos( 



— f°, dcos9- 



dBr 



Pk \ f s F \ 
K~P1 



cos 5, 



KP 
st 



(2m* | A p\ | l|^ P | 2 t 
\ Zs 3 J m? ~T 2 I Ff p I 



(78) 



pKP 

with 5£ — arg.(-p^p-). As we can see from Eq. (1781) . the forward-backward asymmetry is 
determined by the ratio of the scalar and the vector form factors. It is also proportional to 
cosine of the strong phase shift S^ p . The forward-backward asymmetries for Kir and Kr\ 
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cases are shown in Fig. [TTJ As can be seen in Fig. [TTJ the forward-backward asymmetry for 
Kir case is large below K* resonance and reaches to 70%. Here the decay distribution for 
t~ — > K s 7i^u is identical to that of r~ — > K~7r°v by taking the limit for ex of K°K° mixing 
zero. In Fig. [TTJ we have evaluated the forward-backward asymmetry for r _ — > K~tt°i> as 
that for Kit case. 




-1,0 

600 800 1000 1200 1400 1600 1800 
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FIG. 11: The predictions of the forward-backward asymmetries of r — > Kins and r — > Kr}v decays. 
The solid line and dashed line correspond to the forward-backward asymmetry of r — > Kttis decay 
and that of the r — > K-qv decay, respectively. 



V. TWO HIGGS DOUBLET MODEL WITH CP VIOLATION 

As an example of new physics beyond the SM, we investigate a two Higgs doublet model 
with explicit CP violation. (See for example, Ref. [39j for a recent review of two Higgs 
doublet model.) The model is classified as type II two Higgs doublet model. Z 2 parity 
is assigned so that only a Higgs doublet $ 2 is coupled to up type quarks and another 
Higgs doublet $i is coupled with down type quarks. For the charged leptons, they have 
Yukawa couplings with the same Higgs doublet which the down type quarks interact with. 
Z2 symmetry is softly broken in Higgs sector. By taking the soft breaking mass squared 
parameter small, one can naturally obtain the large ratio of vacuum expectation values 



(VEVs) of two Higgs doublets . The idea of Ref. 40J is that this large ratio of the Higgs 
VEVs is the origin of the isospin breaking of the third generation of the quarks. In such 
model, the Higgs with small VEV has enhanced Yukawa couplings to down type quarks and 
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charged leptons. Therefore, r lepton and bottom quark can be good probes investigating 
the extra Higgs doublet with the small VEV. 

The well known effect of CP violation of the two Higgs doublet model is CP even and 



CP odd Higgs mixing 
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42j . In the large limit of the ratio of Higgs VEVs, among three 
neutral Higgs, the SM like CP even Higgs is decoupled from the other two Higgs bosons. 
Therefore, in good approximation, CP even and CP odd Higgs mixing occurs among two 
Higgs bosons in the sector of the Higgs with the small VEV. We investigate how the CP 
violating mixing of the neutral Higgs sector leads to some observable effect on charged 
Higgs Yukawa coupling. We also explicitly show how it generates the direct CP violation 
of t decays. For this purpose, we compute one loop corrections to masses of the charged 
leptons and down type quarks. One finds the one loop corrected mass is flavor diagonal and 
a small CP violating chiral phase due to the CP even and CP odd Higgs mixing is generated. 
To remove the phase of one loop corrected mass, one needs to carry out the chiral rotation. 
After the chiral rotation, CP violating phase in charged Higgs sector arises. The phase is due 
to the CP violation of Higgs sector which is the different origin from Kobayashi Maskawa 
phase {43 1 . 

After all, the relative CP violating phase difference between the charged current interac- 
tion of W boson and charged Higgs interaction arises as, 

C ~ Vli^W^ + VZT R H + e~ 2 ^. (79) 

The phase <p T vanishes if CP even and CP odd Higgs mixing angle 9 ah vanishes. The 
phase <f) T can be measured by direct CP violation of decays. The decays go through 
the intermediate states W~ and H~ which are converted to a common hadronic final state 
(K,ir). Schematically, the process goes as, 

f v L + W~* ] 

t->< } ^ R-jv + v. (80) 

( u L + H-* J 

To measure the phase <j) T , the angular analysis of the decay distributions of r — > Kirv is 
useful. The direct CP violation arises in the interference of two amplitudes with both weak 
phase difference and strong phase difference. In the r — > Kttu decays, the interference of 
two amplitudes with different angular momentum of K~it° , i.e., I = 1 and I = can take 
place. The difference of the angular distribution of and its CP conjugate 
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t + — > K + 7T°i? is sensitive to the CP violating phase described above. As we have shown in 
[nj ], the forward-backward CP asymmetry is a good observable for the CP violation. 

The Higgs potential of two Higgs doublet model with softly broken Z 2 symmetry is given 

as, 



tree 



J2 + y ($!^) 2 ) - m!(${$ 2 + h.c.) + A 3 ($ t 1 $ 1 )($ 2 ,$ 2 ) 

i=l,2 ^ ' 

+A 4 |$ t 1 $ 2 | 2 + ^As^^*^!) 2 + e-^(^ 2 f 



(81) 



where under Z 2 transformation, the Higgs fields transform as, 



$1 



$ 2 -)■ $2- 



(82) 



#5 is a CP violation parameter of Higgs sector. One may write the vacuum expectation 
values with three order parameters 441 ]. 



(*i> 







V% lcos/3 



($2> 



-iff 



V2 I sin/3 



(83) 



The three order parameters are determined by the stationary conditions. For large tan /3, 
the solution can be written approximately as, 



v 2 ~ 



2m 2 



cos /3 



{m 2 + f (A 3 + A 4 )} cos 0' + f A 5 cos(# 5 + 6') 
sin(fl 5 + 0') A 3 + A 4 



sin 6 1 ' 



Ar 



(84) 



where only the leading terms with respect to the expansion of the soft breaking parameter ^§ 
are shown. When 8 5 is not vanishing, the neutral Higgs bosons with definite CP parities, i.e., 
CP even (H) and CP odd Higgses (A) are not mass eigenstates. Their mixing angle is 
sensitive to the CP violation of the Higgs sector. In large tan/3 limit, the mass matrix of 
the three neutral Higgs becomes, 



4u 0^ f h\ 



4 



{h,H,A) 



a 2 2 a 23 
a 23 CJ33 



h 
H 



(85) 
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where ai 2 and are subleading of the expansion of cos/3 and can be neglected in large 
tan (3 limit. Therefore in the limit, one can simply diagonalize 2x2 matrix. For the purpose, 
one introduces the mixing angle 9 ah, 

H = cos 6 ah H 3 + sin 6 AH H 2l 

A = cos 9 ah H 2 - sin 9 ah H 3 , (86) 
where H 2 and if 3 are mass eigen states. The other matrix elements in small cos /3 limit are, 

a u ^ 2A 2 , 

i 

sin 9' 



1 sin(#5 + 9') , 1 , 

a 33 ^ <! V„, ; - cos(0 5 + 29') \ A 5 , 



a 22 - Los9 5 + cos29'^±^)x 5 , 
[ sm 9' J 



a 23 ~ -A 5 sin(0 5 + 20')- (87) 
Then one finds the mixing angle is given by, 

9ah = j + 9'. (88) 

In the same limit, the masses of all the Higgs bosons are; 

M 2 M' 1 XT 2 

£ mass = -?fh* - - - Ml + H+H-, (89) 

with, 



Ml 








M 2 H+ 


= m\ + 


^3 2 
2 ' 






= m\ + 


v 2 

— (A 3 + A 4 — A 5 ), 




M% 3 


= m\ + 


v 2 

y(A 3 + A 4 + A 5 ). 


(90) 



Using the relations in Eq.( |8^|) and the mass formulae of Higgs bosons in Eq. (l90|) . one can 
write the formulae cos (3 and 9 ah as follows; 

cos/3 = m " (91) 

Ml cos 2 9 AH + ML sin 2 9, 

M 2 h^ V, 

-2 
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/ Hz COS- UAH + M- h 2 sm "AH 

I 2 9 " 

9 ah = arctan ( — f^- tan ) . (92) 




FIG. 12: The CP even and odd Higgs mixing angle Bah as a function of CP violating parameter 
#5. The thin solid line, the dashed line, and the thick solid line correspond to = 2, 1.5, and 
1.1, respectively. 



In Fig. 12, we have shown the mixing angle 9 ah as a function of CP violating parameter #5 of 
the Higgs potential as given by Eq. (l92l) . One can see when the mass splitting of if 2 and if 3 
are large, 9 ah tends to deviate from the line of 9 ah = ^f, which leads to 9' is non- vanishing. 
Next we compute the one loop corrected mass due to if 2 and if 3. Yukawa couplings of them 

.'■ \ 

/ s 
/ s 
/ \ 
/ \ 
/ \ 
/ \ 

! I 

FIG. 13: One loop corrections to the self energies of down type quarks and charged leptons due to 
neutral Higgs exchanges. 

to down type quarks and charged leptons can be written as, 

C Y = -1 [tan/^^se-^^m^ + ^e-^^m^O] 
v 

+ El [tan/^e-e-^^m^ + J^ 9 ^ m dl di)} ■ (93) 

Note that the Yukawa couplings of if 2 and H3 have an enhancement factor tan/3. The CP 
violation of the Yukawa couplings are written in terms of the chiral phase, e~ tl5dAH . One 
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defines the one loop corrected masses for down type quarks and charged leptons as, 





f rrii tan (3 ^ 2 " 


I Hoop 





£mass = -lLiMnl Ri - d Li M di d Ri + h.c. (94) 
The corrections are evaluated by computing Feynman diagrams Fig. 13 and the result is, 

-(Zi-l^ + iZ^Zi-^rm, (95) 

where % denote the charged lepton or doun type quark. We have ignored the finite contri- 
bution suppressed by a factor of -rA- and Tjj-- In MSbar scheme, the counter terms are 
determined as, 

Z, - 1 = -(Z mi - 1) = JL (I _ 7 + login) . ( 96) 

Therefore the one-loop corrected masses are finite and are given by, 



4nv J M 2 Hz 



i 



In order to remove the phases of the one loop corrected mass, one need to perform the flavor 
diagonal chiral rotation, 

dm ->• d Ri e- 1 ^ , d Li ->• due* 4 ** , (98) 
where the phases 0, and 0^ are given by, 

sm20 A ^^/J ln^ 
^("W^j cos 2^ In ^ 



tan 20, = ; ; ;"" 7 (99) 



. „ /m d tan/3\ 2 , M% 



tan20 di = ; ^ (100) 

cos 29 AH In M - 



In Fig. 14, we have shown T for different ratios of the Higgs mass j^jf-- The larger ratio 
leads to the larger value of T . 
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FIG. 14: 



as a function of #5. The unit of the angles is in degree. The thin solid line, the 



dashed line, and the thick solid line correspond to -p^ 1 = 2, 1.5, and 1.1, respectively. 



Now we study the effects of the CP violation of the Higgs mixing on r lepton decays. 
The effective four Fermi interactions from the SM contribution and from the charged Higgs 
exchange are given by, 

C cc = 2V2G F V* ; 
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m T mdi tan 2 j3 



101) 



where the relative phase 4> T — of charged current interaction due to W~ exchanged and 
charged Higgs interaction H~ arises. 

The forward-backward CP asymmetry in the two Higgs doublet model can be obtained by 
replacing the SM scalar form factor with the one including the charged Higgs contribution 
in Eq.flTS]), 

,Q 2 tan 2 /3\ KP 



7 KP 

SNew 



-2i(<j> T —(f> 3 ) _ 



Ml, 



(102) 



By comparing the forward-backward asymmetry of r and r + , one obtains the direct CP 
violation fic| . 



A FB (r- ->■ K~Pv) - A fb (t + -> K+PV) 

Q 2 tan 2 $ 



VK 



\F$*\ 



-2sin5f p sin{2((/) r -0 s )}- 



M 



H+ 



(103) 
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where P = 7r°, rj. 




200 400 600 SOO 1000 1200 



M H2 (GeV) 

FIG. 15: The constraints on T parameter in (Mh 2 ,Mh 3 ) plane. The gray shaded regions are 
excluded. We choose M H + = 600GeV. The upper bound on T^ ew is 0.15 and the lower bound is 
—0.068, respectively. The thick solid line, the dashed line, and the thin solid line correspond to 
-jg^ 2 - = 1.1, 1.5, and 2, respectively. 



To predict the CP violation, we take account of the constraints on the mass of the charged 



Higgs, tan/3 and the ratio of neutral Higgs masses jj^-- The lower limit of the chargedJHiggs 
mass is given as M#+ > 295GeV obtained from B — > X s ^ 
B ->• Dtu 



45rj48| . Using B ->• tv 49] and 



50| , the lower limit of the charged Higgs mass is constrained as Mh+ > 500GeV 



for tan/3 



40 48 . The ratio j-^- of the neutral Higgs masses can be constrained from T 



parameter. T parameter of the present model is computed as 511 ]. 

1 



New 



16ttM^s 2 w 



F(M H+ , Mjj 3 ) + F(M H + , Mjj 2 ) — F{M H3 , M H2 ) 



(104) 



where F(m a , m h ) is given by, 



F(m a ,m b ) 



(105) 



In Eq.QH), we take the limit ; Z 3 ->■ f • From Eq.(10.61) of Ref. [521], T New = 0.03 ± 0.11 
for the SM Higgs boson mass Mh = 117GeV case. We shift the SM reference point for the 



Hiegs mass to M h = 126GeV [53j, which amounts to the shift of T New is log iff ~ 0.01. 



Therefore we adopt the following value for TN ew , 



T New = 0.04 ±0.11. 



(106) 
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With M H + = 600GeV, the constraints on (M# 2 , Mh 3 ) plane are shown in Fig. 15. 

In Fig. 16, the forward-backward CP asymmetry in Eq. (I103[) is shown. We neglect (f> s in the 
numerical calculation. We choose the charged Higgs mass M H + = 600GeV and tan/3 = 40, 



and 



M 



»2 



1.1, 1.5, and 2 which satisfy the constraints studied. The CP asymmetry is as 



small as 10 6 ~ 10 7 . Comparing the present result with the one with the two Higgs doublet 



model without natural flavor conservation [lO 
model because it is the one- loop effect. 

4x10 6 



the asymmetry is much smaller in the present 
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FIG. 16: forward-backward CP asymmetry Afb — Afb a s a function of hadronic invariant 
mass \J Q 2 . The solid lines correspond to the ones of r — > Kttv decay and the long dashed lines 



correspond to the ones of r — > K-qv decay. The gray lines correspond to 



lines correspond to 



0.02° and the black 



-0.02°. We choose the ratio of the neutral Higgs mass as = 2, and 



assume the charged Higgs mass as M H + = 600GeV. 



VI. SUMMARY AND CONCLUSION 

Now we summarize our results. For the form factors calculation, we apply the chiral 
Lagrangian including the vector resonance for the computation of the scalar and vector 
form factors of r — > Ki^{j])v decays. 

• We present new counterterms related to the vector mesons and 7/0 in one-loop level of 
pseudoscalar mesons and show how one can perform the renormalization in a system- 
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atic way for the diagram with arbitrary number of the loop. 

• By using the propagator with the one loop corrected self-energy of the vector mesons, 
one can reproduce the vector meson intermediate states. 

• We fit our theoretical curve of the hadronic invariant mass distribution with that 
obtained by Belle. By tuning the parameters, we demonstrate that one can fit the 
hadronic mass distribution up to ~ 1300 MeV well. Between 1300 MeV and 1500 
MeV , our prediction is slightly lower than the experimental data. We also compute 
the branching fraction r — > Kr/v, which is consistent with the experimental one. 

About the CP violation of the two Higgs doublet model, we study the CP violation of the 
Higgs sector. The model was invented to explain the large isospin breaking of bottom and 
top due to large tan (3 ~ 40 [40| . CP violation of Higgs potential leads to the mixings of CP 
even and CP odd Higgs. The Yukawa couplings of down type quarks and charged leptons 
with the neutral Higgs of the second Higgs doublet are large and are CP violating. We found 
that; 

• The CP violation in the neutral Higgs sector leads to the CP violating effect on the 
quarks and leptons mass matrices through one loop corrections. 

• After removing the CP violating phases in the mass matrices, one obtains CP violating 
phases of the charged Higgs couplings to the down type quarks and charged leptons. 

• The effect is studied in the forward-backward CP asymmetry of r — > Kirv decay. The 
order of the asymmetry is 10~ 6 ~ 10~ 7 . The smallness of the asymmetry comes from 
the fact that the CP violation is loop induced effect. 
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Appendix A: Derivation of the counterterms using background field method 



We give the outline of the derivation of the counterterms. To derive the counterterms, 



nc 



d method so that the calculation of the counterterms is consistent 



we use the background 

with chiral symmetry 34] . 351 ] . For the purpose, we first write the chiral Lagrangian in 
terms of the fields which are decomposed into the background fields and the quantum fields 
based on Eq.flTJ). We decompose the fields into the background field and quantum fluctuation 
as follows; 

£l = £exp(zy), £ R = ^exp(-zy), 
fjo = f] + A . 



(Al) 



where £ = exp(iy) and 7r denotes the background pseudoscalar octet fields. A denotes 
the quantum fluctuation. r/ is the background field for singlet pseudoscalar and A is its 
quantum part. We also introduce a± and a defined as 



a, 



2i 
2i 



(A2) 



Using the notations given above, we write the Lagrangian including the background field 
and the quantum parts, 



M 2 

I = / 2 Tr(a ± X) + BTri&Mh + &M£ R )] - ^ 2p Tr(^M£ L - &M£ R ) ■ f)a ~ -fvl 



+M, 2 7 Tr 



a 



M\2 



(A3) 



If we suppress the quantum fluctuation as A — > and A — > 0, then Eq. flA3p equals to 
Eq. (TTT) . We note that under the chiral transformation, £ transforms non-linearly as, 



e = gUrf = h& 



while A transforms linearly as, 



A' = hArf. 



(A4) 



(A5) 
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We also note a and a± transform as, 

1 



-/ 
a 



' M = ha^lv + -hd,Jv , 



a 



i 

-i 



ha^rf. (A6) 



We treat the vector meson V as the background field and it transforms as 

VI = hVfl + -.hd^. (A7) 
Given the transformations in Eq. ( 1A4I) -Eq. ( 1A7I) . the Lagrangian of Eq.( 1A3l) with the chi- 



ral breaking term M replaced with the spurion fields [34[ is invariant under chiral 
SU(3)l x SU(3)r transformation. We next compute the one-loop corrections and iden- 
tify the divergence. For the purpose, one expands £ up to the second order of A and 
A . 

/ M 2 \ 

C = £ + TtD.AD^A + M - -^J Tr[A, a Ul ][A, 
+ ^Tr(^-^)[A,/^A] 

- ^Tr( X+ A 2 ) + ^Tr^A 2 )^ + 2^Tr( X+ A) A 

+ -d^A°d^A° - M A ° 2 , (A8) 



where \± is defined as, 



X- = i(£M£-?M?). (A9) 



Since the background fields (£, r/ ) satisfy the equations of motion 



M 2 ol^ 1 

D"a ± „ + l j^[ vfl ~ —, a±„] = p( B X- + 92 P X+Vo 

□ + M 2 )r/ + g 2p Tr(X-) = 0, (A10) 



the first variation with respect to A and A vanishes. Introducing the octet component 
field A a as A = X^o=i A a T a , one can write the quadratic part of the action in terms of the 
quantum parts A A = (A , A a ) as, 

S = S-l [ d 4 xA A D AB A B , (All) 



2 
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D AB is a diffrential operator and a 9 x 9 matrix for the nonet space. 

D AB = A D + M o) 

I M a0 D ab 

lor" 6 / \ o ^ 



with, 



where, 



and, 



We also define, 



and 



D ab = {d^) ab + a ab - (v 2 ) a \ 



(d^A) a ^(d,A) a + vfA b , 
(d,Ay = d,A a + TfA b , 

T ab = _2iTr[T a ,T 6 ]aM> 
V ab = _2iTr[T a ,T 6 ]\/ M , 



Ml , 1 

v 



2gf 2 v g 



a ab = ^Tr( X+ T a T b ) - ^Tr( X _T a T% 



-2 1- 



g 2 Pj 



M a0 = M 0a = _^P Tr ( x+T a) 

The effective action including one loop corrections is given by 

S eff = S + AS, 

AS = -TrLnD AB . 
2 



By introducing, 9x9 matrix, 



a 



w i Ml M° b 

M a0 a ab 



40 



one can write D AB as; 



D 



AB 



{d^) AB + a 



AB 



,2\AB 



(A20) 



The diverge nt p art of one loop correction can be easily computed with the heat kernel 



method 



341 ] [361 ] . The counterterms can be also obtained with, 
S c = -AS\dw = ^ J d 4 xTi[a 2 (x)], 



where a^(x) is given by, 



c,,U) = + Ud^d v ][d^d v ] + \[d„ 



(A21) 



(A22) 



2 12 L JL 6' 
where a = a — v 2 . The trace for 9x9 matrix can be converted to trace for 3x3 matrix 
which leads to the conterterms of Eq. f[2"Tj) . 



Appendix B: 1 loop functions 

Here we summarize the one loop functions which appear in Eq. (1301) . 

d d k 1 



(2ir) d i k 2 — m 



2 ■ 

P 



qp _ I _£k (Q - 2k), 

X V " / (0^\dA S(U _ n\2 _ m 2 , \(1.2 _ ^2 V 



(2ir) d i {(k - Q) 2 - m 2 Q }(k 2 - m 2 P ) 
QP f d d k (2k-QM2k-QW 

" ~ J {2-n) d i Uk- Q) 2 - m 2 \{k 2 - m 2 P ) ' 1 ' 



i d 

(2n) d i {(k - Q) 2 - m 2 Q }(k 2 - m 2 P ) 
By carrying out the loop integral, one obtains 



I p = —2m 2 p \ — 2f 2 /i P , 



+ y-Q, ( - 2AS P0 - 2(^ Q + » P )f 2 ) , (B2) 



Jpq is defined in Eq. flClOj) and Eq. (1C12|) . Hpq is defined in Eq. fl33|) . By taking account of 



T] and r/s mixing, one defines I m and XJf K (X, = x^, J^) as, 

I m = I v cos 2 6» 08 + I v , sin 2 6 08 , 

XVS K = X «K cos 2 ^ + x rfK gin 2 ^ (fi3) 
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Appendix C: Two point function of vector mesons 



Let us determine the coefficient of the counterterms C\,C2 and Zy from the renormal- 
ization for self-energy of vector mesons. The vector mesons couplings with pseudo scalar 
mesons in C are, 



C 



VPP 



_Ml 

gf 2r i 

Myl 

+ 2p A 



Tr(iyA,^A]) 



?T d TT 



V2 



(CI) 



The quantum field for pseudoscalar octet A is denoted by n. One can parameterize the 
inverse propagators of vector fields as, 



A y f + B V Q»Q U . (C2) 
We study the two point functions for p + and K* + mesons. 

A v = M v -Q 2 + 5Ay, 

By = l + 5By, (C3) 

where 5 Ay and 5B V denote the one loop corrections including the contribution from the 
counterterms. For the p meson, they are given as, 

6B p = Z\{p) + {^j (Ml + l -M^ , 

5A P = -^p-(^) 2 (^ + ^x) f 2 

+C r Mml + C r M(2m 2 K + ml), (C4) 

(P = 7r, K) are the loop functions for ir mesons and K mesons 



where/i P = 3^1n^. M r p 
defined as, 

1 

12 



M r P = — 

r 1 O 



1 - 



4m^ _ 1 j mp 
IQ-ir 2 p? 



Q 2 



48tt 2 



4m 2 



Jp = < 




l-l/l- 



4m 



4m 2 „ -. 

-t^j- — 1 arctan 



(C5) 



4m p_i 
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Zy(/i),C[([i) and C£(/i) are finite parts of the renormalization constants defined by, 

r2 \ 2 



Zy = Zy(fi) - 

Ci = c 'M - 
c 2 = c r M - 
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1287T 2 


\gP 


3 ( 


(Ml 


1287T 2 1 


KgP 


1 ( 


'My 


1287T 2 » 


\9P 



(C uv + 1 - In// 2 



(Cuy + l-lnfl 2 ), 



(C uv + 1 - In// 2 



(C6) 



with Cuv is the divergent part of the dimensional regularization, 



Cuv = - - 7 + ln(4vr) 

e 



(C7) 



where e = 2 — | and 7 is Euler constant. The self energy corrections to K* meson are given 
as . 



8B K * = Zy(») + \ 



2 \ 2 



5A* K 



-Q SB K * + - 



3 / M v x 2 



4 \gP 

+C r Mm 2 K + C r M(2m 2 K + ml 



(C8) 



where Mpg and Lpq are the same functions as the ones defined in Ref. 34] 

I / S 



i2g 2 

kpQ 



SQ 4 



-, ■ , -po m 2 p m 2 n m 2 n 
32tt 2 \ ApQ Ap Q m 2 P 



6 288vr 2 ' 



A 2 
4s 



(C9) 



masses m P and rriq, Above the threshold; Q 2 > {m P + mq) 2 , it is given by, 



where kpQ = — a^ -17 ■ </pq is a one loop scalar function of pseudo scalar mesons with 



Jpq{Q 2 
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32tt 2 



2 + 



A, 



In 



Q 2 m 2 p A PQ m 2 p 

2 i \2 A2 



z/pQ (Q 2 + v PQ ) 2 ~ Ap Q 
Q 2 n (Q 2 - v P q) 2 - A 2 PQ 
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\ VPQ 

i6vr g 2 ' 



(CIO) 



where, 



j pq — Q 4 ~ ^Q 2 ^pq + ApQ> 



(Cll) 
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while below the threshold (mp 

1 



m Q ) 2 < Q 2 < (m P + m Q ) 2 , 



Jpq(Q z 



32tt 2 
- 2 
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arctan 
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+ arctan ■ 
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PQ 



with E = m 2 P + rrig and Apg 



mr> — m 



Q- 



(C12) 



Appendix D: Self-energy corrections for rj — rf sector 



Here we have expanded the interaction terms which are needed to compute Zs8 ; ^Mq 8 anc ^ 



<5Mf 8 in Eq.(j6]) as well as the wavefunction renomalization constants Zr and z m in Eq.(|3T 
within one loop approximation. The relevant part of the Lagrangian for the calculation is 



C 



1 - 
IB 



M 



9 2 P 



P 



2 Tr( X+ A 2 ) + ^Tr( X _ A 2 )^ + 2^Tr( X+ A)A° 
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^Tr(A5 M 7rA^7T - A 2 c>^7r) - glTr({7r, X }A 2 )r/ 

1 



92p 



Tr({7r, {tt, X }}A) A - — -Tr({7r, {vr, X }}A 



;di) 



In terms of the component fields, the Lagrangian is, 



-^drj 8 drj 8 (k°k° + K + k-) - V8Vo ^L x 



5m 2 — %m\r 3m 2 — 4m 

K VsVs + - 



_ K (k°K + K+K-) + ^(7r°7r° + 27T+7T- 
2\/3 2^ 
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-VoVsiVs 1 



2Bf 
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lQm 2 w — 7m 2 ~ Sm 2 ^ — 3m 2 

— vl + 



* (k°K° + K+K-) + ^(7r°7r° + 27T+7T- 
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• w , 3. 2 K°K° + 2K+K 
+ -dK + dK~(-- V 2 8 
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7T 7T 



2 4 



(D2) 



We denote the quantum field for pseudoscalar octet A as 7r and singlet Ao as ?fo. i]o,K,7r and 
7/g denote the background field. The relevant counterterms can be extracted from Eq. fpJT]) 
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as, 



8L 5 4m 2 K — m 2 _ 8L 6 2 2 .4m|- — m 2 2 



16L 8 (H-mJf, 2 

/2 1 6 + 3 

4T 3 g 2p , 2 , 2n , &T 5 g 2p 2 , . 

Tffi/ ^ x + m w )r] r] 8 + —-^=m K A K7r r] r] 8 . (D3) 

Appendix E: The form factors for r — >■ fTryf decay and r — >■ Ifr/V decay 

In this appendix, we give the equations of the form factors for r — > Kr/v and r — >■ Krjv. 
In the following equations, 5 A and SB imply SAk* and 5Bk* , respectively. In this appendix, 
we give the equations of the form factors for r — > Krju and r — > Krj'v. The vector and scalar 
form factors are given as the sums of the contribution of 1 PI diagram and K* resonance 
contribution. 

W Kr > — T? IPI _l_ T? K * 
r V ~ r VK-q "I" r VK v i 

= Fsk v + F s % (El) 

The contribution of the 1 PI diagrams for r — > Kr/u form factors is computed as, 
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(K + r]\u L -/ fl s L \0}\ 1PI = cos9 08 -j= x 
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Using Eq.( 1E2|) . the form factors of 1 PI part is derived as, 
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The decay amplitude of the process K* — >■ is given as 
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with 



E Kr) = V3 cos 6*08 
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Using 



v 5m x ~ 3m* / A K7r - 1 - 2 A w - 2 
32^ Q 5 " 3 (J ^' COS ^ 08 + Sm ^ 08 

the .fT* — > i^r? decay amplitude, the contribution to the form factor is given by, 



SKl > " Q 2 M V + 5A + Q 2 5B' 

The form factor for r — >■ X77V is also given as, 
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The decay amplitude of the process K* — >■ K?/ is given as, 



T M (ir + ->■ #+7/) = Ekti'Ill + Q„A Kv ,F Kv ,, 
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with 
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Using the K* — > Krj' decay amplitude, the contribution to the form factor is given by, 
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